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Abstract 

In this paper, we consider the boundedness of solutions for a class of impact oscillators 
with time dependent polynomial potentials, 

In 



x + x 2n+1 + ^ Pi {t)x l = 0, for x(t) > 0, 



i=0 

x{t) > 0, 

H t o) = ~ i ( t o)> if x(t ) = 0, 

where n G N+, pi(t + 1) = pi(t) and Pl (t) e C 5 (M/Z). 

Keywords: Impactor oscillators; Boundedness of solutions; Canonical transformation; Time- 
dependent polynomial potentials; Moser's small twist theorem. 



1 Introduction and main result 

In Ref. p], Dieckerhoff and Zehnder proved the boundedness of solutions for a time-dependent 
nonlinear differential equation: 

In 

x + x 2n+1 + ^2p i (t)x i = 0, n>l, (1.1) 

i=0 

with pi{t + 1) = Pi{t) and Pi(t) E C°° and asked whether or not the boundedness phenomenon 
is related to the smoothness of Pi{t). Then Laederich and Levi [2] relaxed the smoothness re- 
quirement of Pi(t) to C 5+e with e > 0, Yuan Od] relaxed the smoothness requirement of Pi(t). 
to C 2 , etc. 
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Motivated by [TJ [21 0] , we are going to consider the problem of the boundedness for Eq. (jl.ip 
with impacts: 

( 2n 

x + x 2n+1 + Y J Pi{t)% i = 0, for x(t)>0, 

•=o (1.2) 

s(t) > 0, 

> x(4 ) = -x(to), if x(i ) = 0, 

where n G N+, p;(i + 1) = and pi(t) G C 5 (R/Z). 

The system is included in the system of impact oscillators given by 

x + V x (x, t) = 0, for x(t) > 0, 

x{t) > 0, (1.3) 
Hta) = -»(*o )> if x (*o) = 0, 

which serve as models of dynamical systems with discontinuities [5]. From the viewpoint of 
mechanics, this equation models the motion of particle attached to a nonlinear spring and 
bouncing elastically against the fixed barrier. The system of this form also relate to the research 
of the Fermi accelerator [6] , dual billiards [7] and celestial mechanics [8] . 

The nonsmoothness caused by the impact limit the applications of many powerful mathe- 
matical tools. However there are also many interesting papers on the impact oscillators, see 
[9]-|16] and their references. 

To overcome the problem of nonsmoothness, Qian and Torres [9] and Qian and Sun [10] used 
success mapping in their papers and proved the existence of invariant tori by a variant version of 
Moser's small twist theorem by Ortega [17] . In Zharnitsky [11] , Wang [121 113] . they exchanged 
the position of angle and time variables of the Hamiltonian function which only C° in angle 
variable, and finally obtained the boundedness of solutions by Moser's small twist theorem. 

In [12], Z.Wang and Y.Wang studied the following impact oscillator: 

x + x 2n+1 = p(t), for x(t) > 0, 

x(t) > 0, (1.4) 
^o") = -^o )' if x ( t o) = °' 

where n G N + , p(t) is a C 5 periodic function with 1, and the term x 2n+1 models a hard spring. 
They proved the Lagrangian stability for system (II. 4p by exchanging the roles of angle and time 
variables and using Moser's small twist theorem. 

In this paper we will extend (|1.4p to the case of time-periodic polynomial potentials with 
constant leading coefficient as (jl.2p and use similar method to prove the boundedness of solutions 
for (jl.2p . More exactly, we obtain the following conclusion: 

Theorem 1 Every solution of (1.2) is bounded, i.e., x(t) exits for t G K and 

sup(|x(t)| + \x(t)\) < +oo. 
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The idea of proving the boundedness of solutions of (jl,2p is as follows. First of all, we 
describe (|1 .2|) by a Hamiltonian function H^(p,(j),t) (see (|2.6p ) in action-angle variables defined 
on the whole space R + x S 1 x S 1 by means of transformation theory. Due to the existence 
of impact, H$(p,(j>,t) is only continuous in (j), but is C 5 smooth in p and t. Exchanging the 
roles of the variables <f> and t outside of a large disk D r = { (p, <p) , p < r} in (p, </>)-plane, (2.9) is 
transformed into a perturbation of an integrable Hamiltonian system H±{1, 9, r) (see (3.3)) which 
is sufficiently smooth in / and 9. The Poincare mapping with respect to the new time r is closed 
to a so-called twist mapping in (R + x S 1 ) \D r , and satisfies Moser's invariant curve theorem 
after a scaling transformation. Then Moser's theorem guarantees the existence of arbitrarily 
large invariant curves diffeomorphic to v = const, over v = in (v, #)-plane. Go back to the 
equivalent system (2.1), every such curve is a base of a time-periodic and flow-invariant cylinder 
in the extended phase space (x, y, t) E M + xIxK, which confines the solutions in the interior and 
which leads to a bound of these solutions if the uniqueness of initial value problem holds. Note 
that, from the proof below, system (2.1) is equivalent to a smooth system (2.9) 
which has uniqueness. Hence system (2.1) has uniqueness. 

Remark 1.1 The time- dependence of the stiffness coefficient Pi(t) of the spring can be produced 
by periodic changes of the temperature or other physical variables. Following Laederich and Levi 
and Z. Wang and Y. Wang [Wj , we assume Pi(t) £ C 5 . 

Remark 1.2 Compared with (1-4), the potentials of (1.1) become more complicated which will 
make more difficulties in the proof of Theorem 1. 

Remark 1.3 The main difficulty in this paper is also the nonsmoothness caused by the impact. 
Similar to UStf . we will also exchange the roles of angle and time variables of the Hamiltonian 
function H-i(p,4> : t) (see below), and then obtain the boundedenes of solutions by Moser's small 
twist theorem. 

The rest of this paper is organized as follows. In section 2, we will make some action-angle 
transformations to transform the system into an equivalent Hamiltonian H^(I, 9,t) (see below) 
which is defined in the whole plane IR + x S 1 x S 1 and C 5 in /, 9, but only C° in r. In section 
3, we do more canonical transformations such that the Poincare mapping of the new system is 
closed to the twist mapping and then use Moser's theorem to complete the proof of Theorem 1. 
Finally, we will give the proofs of some technical lemmas in section 4 and draw a conclusion in 
section 5. 

Throughout this paper, we denote A < .B if A < C\B, and A > .B if A > C2B , where 
C\,C2 are positive constants. 

2 Action-angle variables 

In this section, we will make some action-angle transformations to transform the system into 
an equivalent Hamiltonian H^(I, 9, r) (see below) which is defined in the whole plane R + xS 1 xS 1 
and C 5 in /, 9, but only C° in r. 
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Without impact, Eq. (jl.2p is just Eq. (|l,ip which is of second order and equivalent to the one 
order system 

x = y, 

y = -^ +1 -f>(t)*\ (2J) 
(2.1) is a system defined in the whole phase plane XOY, which has a Hamiltonian function 

H(x, y, t) = -y 2 + — L_ x 2n+2 + V ElM^ (2.2) 
v y ' 2 y 2n + 2 ^ i + 1 v ; 

with the symplectic form dx Ady. 

In order to make an action-angle transformation, we consider the system 

x = y, 

y = - x 2n+1 . 

i.,,2 _l_ 1 x 2n+2 



(2.3) 



with the Hamiltonian Ho(x,y,t) = -^U + 2n +2 

Clearly, Hq > on R 2 except at the only equilibrium point (x, y) = (0, 0) where Hq = 0. 

Let (C(t),S(t)) is the solution of Eq.(2. 3) satisfying the initial condition: (C(0), 5(0)) = (1,0). 
Let To > be its minimal period. 

From Eq.(2.3), we can find that C{t) and S(t) satisfy: 

(1) C{t), S{t) € C°°(M), and C(t + T ) = C(t), S(t + T ) = S(t) with C(0) = 1, 5(0) = 0; 

(2) C{t) = S(t) and S(t) = -C{t) 2n+1 ; 

(3) (n + l)S(t) 2 + C(t) 2n+l = l; 

(4) C(-t) = C(t) and S(-t) = -S{t). 

Now we can make coordinate transformation by the mapping \ : M + x 5 1 — > M 2 /{0}, where 
(x,y) = ^i(X,'d) defined by the formula 

x = (a\) a C(VT ), y = (aA) Q 5(tfT ), 

where a = ^r^, (3 = 1 — a, a = ^ are constants. We can claim that is a symplectic 
diffeomorphism from R + x 5 1 onto R 2 /{0}. Indeed, the Jacobian determinant of is 1, so 
is measure preserving. Moreover, since (C(t),S(t)) is a solution of Eq. (|2.3p and has To as its 
minimal period, which concludes that *$>i is one to one and onto. This proves the claim. 
Under Hamiltonian H is transformed into 

BlMt)= ^ + f : ^y^y^f H i ( , 4) 

2n + 2 f— ' i + l 

which is C 5 in A, t and C° in "9. 
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Considering the impact case, the phase space is only a half plane (x > 0) of the original 
phase plane XOY. Since impact, the smoothness with variable $ of H\ is disappeared. We shall 
study another Hamiltonian H%(p,(j),t) as defined below. Let C(i?To), we will find two solutions 
$ = #2 in interval [0, 1). By the symmetry of Hq, we can easily find i?i = 0, $2 = 5 ■ 

Under (A, 0(mod 1)) is mapped into the positive Y-axis, (A, ^(mod 1)) is mapped into 
the negative Y-axis. 

Define a sympletic transformation ^2 : ((fr^p) ^ A), given by $ = ^, A = 2p. And H\ is 
transformed into 

H 2 (p,<j),t) = dp p + y — — , (2.5) 

j-i 1 + 1 

where d=|gj. 

To get a system which is equivalent to (1.2), we could define a new Hamiltonian function by 

to 1 + 1 

where [<t>\ denotes the largest integer less than or equal to eft. And the corresponding Hamiltonian 
system is 

p=^_ : when <j)(t) € (k,k + l),k 6Z, 

^ = -^» when 0(t) € (fc,fc + l),A;GZ, (2.7) 
0(to) = p(*o) = 1™ p(t), when lim (j>(t) = k,k € Z. 

Obviously, #3 is periodic in </> with 1, and C 5 in p, <^> when 4> ^ Z, but only continuous in 
when € Z. In fact, #3 is right continuous when € Z. 

Remark 2.1 T/ie systems mentioned above are equivalent, see \12j and U3\j , we omit the proof. 



3 The proof of Theorem 1 

Now we are concerned with the Hamiltonian system (|2.7p with Hamiltonian function #3 given 
by (12.6|) . To cope with the nonsmoothness in <j) and use Moser's small twist theorem to prove 
the Lagrangian stability, we will exchange the positions of variables {p,4>) and (H^,t) below. 
This trick has been used in [TH1 f!9l [20] . 
From (|2.6|) . we have that 

lim H3 = 00 

p— i-OO 

and 

9Hs _ 0AR 2 -i , ^ «(i + l)^(O(2a)^ +1 )^ +1 )- 1 {g(^(0-[0])^o)r +1 > n 

dp ~ app + ^ i + i >u ' 

' i=0 
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if p > 1. By the implicit theorem, we know that there is a function R = R(H 3 , t, eft) such that 



p(H 3 ,t, 0) = t r*PH> p - R(H 3 , t, </>). (3.1) 



Then we make another transformation ^3 : (p, <j), t) 1— > (J, 0, r) given by 

I = H 3 (p,cf>,t), 9 = t, t = 4>- (3.2) 
This transformation also leads to a Hamiltonian system with new Hamiltonian function 

H 4 (I,,0,t) =d~^I^ - R(I,0,t), (3.3) 

which is C 5 in 6, C°° in / and C° in r. 
And the system is given by 

= ^, r(jtZ, (3.4) 
9{k) = lim 0(r), I(k) = lim I(r), fe G Z. 

T— >fc T— >fc 

Then we have 

Lemma 3.1 For / /an?e enough, the following estimate about R hold ture 

\D{D j e R\ < J5 -1 , for i + /c<5. (3.5) 
To be convenient for readers, the proof is given in section 4. 

In order to use the Moser's small twist theorem, we introduce a new variable v and a small 

positive e by the formula I = j, where v S [1,2], e > and e — > 0, when I — > 00. Then the 
system (|3.4p is equivalent to the following system 

9 = ^d-^e^Wv^' 1 -s^, t(£Z, (3.6) 
9(k) = lim 0(r), = lim v(t), keZ. 

Integrate the system (|3.6p from t = to r = 1, we obtain the Poincare mapping P of the 
form 

vi = v + fi(v ,e ), 0i = e + ^d-^e 1 -^v^~ 1 + f 2 (v ,e ), (3.7) 

where fi(v , 6 Q ) = f Q eff dr, f 2 (v , 9 ) = - f Q e^dr, which satisfy: 
Lemma 3.2 If e small enough, then 

\Di D k e h(vM <J, for j + fe<4; (3.8) 

and 

\Di D k e J 2 (vo,9o)\<J, for j + k < 4. (3.9) 
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Proof. Note / = v/e, we have D 3 v R(v/£,9,t) = e ■ J \Dji?(J, 6, r), and e — > 0, when / — ^ oo. 
Thus, with (|3.5p - (|3.7p . it is readily to get the results, see pQ. We omit it. 

Lemma 3.3 The mapping P possesses the intersection property in the annulus [1,2] x S , i.e., 
ifT is an embedded circle in [1, 2] x S 1 homotopic to a circle v = const., then P(T) n T ^ 0. 

The proof can be found in pQ. 
Proof of Theorem 1. Until now, we have verified that the mapping P satisfies all the condi- 
tions of Moser's small twist theorem HH [22]. Hence for any oj £ (^oT^e 1-5 ?, ±d~v>2 1 ~w) 
satisfying 

> .e|g|~ 5/2 , (3.10) 

there is an invariant curve r of P surrounding vq = 1 if e is sufficiently small. This means that 
the Poincare mapping of the system fj3.4j) has an invariant curve V surrounding I = - and in the 
annulus ^ < I < ~ . Then transform the system (|3.4p back to the system (|2.7|) , we can find that 
there exit arbitrary large invariant tori in (p, <f>, t) space, and this implies that all the solutions 
are bounded. Finally, according to the equivalent of the system and (12. 7)1 . we can imply 
the boundedness of all the solutions of Eq. (|l.ip . The proof of Theorem 1 is completed. 



4 The proof of Lemma 3.1 



Proof. From (13.1)1 and <\3.2h . we have 

p = d~^I^-R (4.1) 



i i 



where a = /3 = 1 - a = ^±±, a = ^ and d = gj. 
When j + fc = 0, substitute p in ()2.6|) . we have 

, i i v5fl ^p i (0)(2a) a ( i+1 )(^^/^ -i?) a ( i+1 ){C(i(r- [r])To)} l+1 , 
/ = d{d--pIWi - RfP + ^ ™l2± — , (4.2) 

and 

^ Pi {e){2a)^ l+1 \d~^I^ -R) a ^{C{Ur- [r])T )} i+1 , _i jl nM , , 
^ — = d(d wlw -R) 2p . (4.3) 

i=o 21 + 

Using the Taylor expansion of / CT (cc) = (1 + x)°" = 1 + ax + x 2 f"(sx)(l — s)ds, where 
f"(sx) = o~(a — 1)(1 + sx) a ~ 2 , and with the constants a, (3, d, we have 

d(d~^I^ - Rf 3 = I(\ - d^rWR)^ = X- 2j3d^I l ~^R + IR U (4.4) 
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where 



and 



•l 

R x = (d^r^R) 2 I fUsd^r^R)(l - s)ds, (4.5) 
o 



1 1 ^2 



^ Pi (9)(2a) a(i+1 \d-^I^ - i?) a ( i+1 ){C(i(r - [r])To)} i+1 
i=Q 

7 - ^ /i 4 (6», t)I% (1 - a(i + l)d^T^R + # 2 ), 

j=0 



where 



where 



2n 



i=0 

1 i i 



(4.6) 



p t (g)(2a)^+ 1 )d-^{C(^(r - [r])^)}^ 1 

Vi{8,V = , (4.7) 

i? 2 = (d^T^R) 2 J f^ i+1) {sd^rWR){l - s )ds. (4.8) 

From (|4T3j) - (|4T8j) . we have 

fl(J, 9, t) = gi(I, 9, t) + 02(1, 9, r) + g 3 (I, 9, r) + g 4 (I, 9, r), (4.9) 



<7i(I,0,r) = £(2/3d*)- Vi(«,r)/^— (4.10) 
i=0 

2n (+i)+l 

<te(I, 0, r) = - Y^(2P)- x ih{0, r)/ £ T i - 1 i?, (4.11) 



2n (+i)+i 

ff3 (I, 0, r) = Y^pdk)- 1 ^ r)I^T- l R 2 , (4.12) 



ff4 (J, 9, t) = (2/3d^)- 1 I^R 1 . (4.13) 
For / large enough, we can prove the following inequality hold true, 

\R\<JV, (4.14) 

from dOt- lOS]) . 

Assume we have proved (|3.5p for j + k < 4 already, i.e., for I large enough, 



\D j I D k e R\ < .n-3, for < j + k < 4. (4.15) 
We are going to prove {D^D^Rl < for j + k = 5. From (Oil , we have 

D^i? = D{D k e9l + Dfzfe + DjD$g 3 + D?Z>^ 4 . (4.16) 
We will give the estimates of every term of the right hand of Eq. ()4.15p as follows. 
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Lemma 4.1 

\D j I De9l\ < for j + k = 5. (4.17) 

Proof. From (j4.10|) . we have 

2n 

Dj^ 5l (/,e,r) = ^(2^*)- 1 J D^/x i (0,r)/ a ^ - 1 , (4.18) 



i=0 

Using Leibnitz's rule, we can easily prove 

g(i+l) 1 a(i+l) 1 



D^i(9,T)I^ x | < .J - ^ iorj + k = 5, (4.19) 



from (14. 7p and the definition of Pi(r). 
Finally, from (|4.16|) . we have 



2n 

i(»+i) 



\D{D k e9l {I,e,r)\ = \Y j m^r 1 D :i I D k e r lt (e,T)I^- 1 \ < .1^ . 



i=0 

Lemma 4.2 



\D 3 I D$g 2 \ < .r^~ j + I~^\D j I D k B R\ 1 for j + k = 5. (4.20) 
Proof. From (|4.1ip . we have 

2n (+l)+i 

DlD k 9 g 2 (I,e,r) = - J^)- 1 !^//^, r)!^^'^, (4.21) 



i=Q 

Using Leibnitz's rule, we have 

»(>-i; I v • , . , ; 



31 + 32 — 3, 
k\ + k 2 — k, 
j + k = 5 

DfD^i^^I^-^DfD^R + ^^I^^D^R. 



31 + 32 — 3, 
kl + k 2 — k, 
j + k < 4 



Then we can prove that 

\D\D k e g 2 \ < .r^~ j + r^\D j jD%R\, for j + = 5, 
from (f4T7|) and P~T5j) . 



(4.22) 
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Lemma 4.3 

< .J5S3~ J + J~as+i|Z)JZ)Ji2|, for j + jfc = 5. (4.23) 



Proof. From (|4.13p . we have 

2n 

7 * — T 1 -i-7 + + l -, 

D^D k 9 g 3 (I,e,r) = ^(2/3d^)- 1 D^ e V(^ r)n (4.24) 

i=0 

Firstly, we give some claims before proving (I4.23p . Similar to |12j Lemma 3.3, we have 
Lemma 4.4 Let F : R — > R and f : M 2 — > R be sufficiently smooth, then for j + k > 1 

DlD k y [Fof] = Y, C sU {D s f F{f)){D{}D k y -f)...{DifD k ^f) (4.25) 

1 < s < j + fe, 
i = Ul, ~,ja), |j| = j! 

fe = .... fc s ), |E| = fe 

where |j| = j a + ... + j s , |fc| = k x + ... + fe a . 
From (14.81). we have 



D j jD k R 2 = DjD k (d^r^R) 2 j f^ i+1) {sd^r^R){l - s)ds. (4.26) 



and 

\^D k ^f^ +l) {sd^r^R){l - s)ds\ < .\D\D^f'^ +l) { S d^r^R)\ 
< .\D j jD k (l + sd^r^R) a - 2 \ = \D j jD k F a (r^R)\ 



Let F a {x) = (1 + sd^x) a ~ 2 , then we have 



D s x F a {x) = -f-F a (x) = C as (l + sd^xy- 2 ~ s (4.28) 
dx s 

where C as are constants. 

_ i 1 

With (j4T25l) and (1X281) and |J af.R| < ±, provide I large enough, we obtain that 
\D\D k [F(r^BR)]\ 

= I E C -J ^W(/-*ll))(2?} 1 i3{ 1 (/-*fl))...(I?}'i3*'(J-*fl))| 



(4.29) 



1 < s < j + fe, 
j = Ul! -lie), iJl = J. 
fe = (fci, fe s ), |fe| = fe 

< • £ |(i^^(r-*ii))|...|(D}'i3}-(i-*fl))| 

1 < S < j + fc, 

j = Oil ■■•7j«)> b' = j> 
£ = (fei, fe s ), |fe| = fe 

For j + k < 4, we have 

{DjD^r^R^ < JS&s-i. (4.30) 
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When j + k = 5, 

|Z^D^(/~^i2)| < J~5S+5~ J + r^+2\DjDlR\. (4.31) 
From (I4.28p -( l4.30p . we can prove that 

|Dj£>£[i^(I-^12)]| 

< j- 2^+2- i + /-2^T2-J + /-2^F5-J + i-ii+2-i + (4.32) 

< J 2«+2 J _|_ / 2n+2 for j + fe = 5. 

Thus, let a = a(i + 1), we have 

\DiDo[Fa(i+i)(I~^R)]\ < J~^~ J ' + r£&\D 3 I D$R\ for j + k = 5. (4.33) 
Using Leibnitz's rule again and with (|4.30p - (|4.33p . we can obtain 

\D j jDlg 3 \ < .r^vi- j + r^\D\D^R\ for j + k < 5. (4.34) 
Finally, we prove that 

2n | C+lJ+1 l 2 

|£^ 53 (/,0,t)| = \Y J (^d^)- 1 D :i I D k e L M (d,T)l' L2 r- 1 R 2 \ < ./^^ + J-J&l [D^Rl 
i=0 

(4.35) 

Lemma 4.5 

|I>fDeS4| < .I^2' j + r^\D j I D 1 lR\ for j + jfc = 5. (4.36) 

The proof is similar to Lemma 4.3, we omit it. 

Now, with lemma 4.1, 4.2, 4.3, 4.5 and (j4.15j) . we have 

\D{D k e R\ 

< Js-J + i~2^+2- J + 7-2^+2 \D 3 I D k e R\ + J2^+2~ J + 7-2^+2 \D\D k R\ + 7^+2^ + J _ 3m5 IDjDjfll 

< .Ja~j for j + fe = 5. 

(4.37) 

Thus, we have \DjDqR\ < ./2~ J provided < j ' + k < 5. The proof of Lemma 3.1 is 
completed. 

5 Conclusion 

In this paper, we study the Lagrangian stability for a class of impact oscillators with time 
dependent polynomial potentials. By exchanging the roles of angle and time, we overcome 
the nonsmoothness on angle variable caused by the existence of impact. Then we obtain an 
equivalent Hamiltonian which is sufficiently smooth in new angle variable by some canonical 
transformations. Thus by direct application of Moser's twist theorem, we obtain the existence 
of invariant curves for the Poincare mapping of the equations, which implies the Lagrangian 
stability. 
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